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Abstract: The public is increasingly concerned about particulate matter pollution caused 
by respirable suspended particles (PMio) and fine particles (PM2.5). In this paper, 
PM10 and PM2.5 concentration are estimated with remote sensing and individual air quality 
indexes of PMi 0 and PM2.5 (IPM10 and IPM2.5) over mainland China in 2010 are calculated. 
We find that China suffered more serious PM2.5 than PM10 pollution in 2010, and they 
presented a spatial differentiation. Consequently, a particulate-based air quality index 
(PAQI) based on a weighting method is proposed to provide a more objective assessment 
of the particulate pollution. The study demonstrates that, in 2010, most of mainland 
China faced a lightly polluted situation in PAQI case; there were three areas obviously 
under moderate pollution (Hubei, Sichuan-Chongqing border region and Ningxia-Inner 
Mongolia border region). Simultaneously, two indicators are calculated with the 
combination of population density gridded data to reveal Chinese population exposure to 
PM2.5. Comparing per capita PM2.5 concentration with population- weighted 
PM2.5 concentration, the former shows that the high-level regions are distributed 
in Guangdong, Shanghai, and Tianjin, while the latter are in Hebei, Chongqing, 
and Shandong. By comparison, the results demonstrate that population-weighted 
PM2.5 concentration is more in line with the actual situation. 
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1. Introduction 

In recent years, inhalable particulate matter has become the primary air pollutant in China [1]. 
Many studies have shown that fine particulate matter concentration has a certain association with 
pollution-related disease mortality [2-5]. Public opinion now pays more and more attention to PM2.5 
because of its connection to health risks [6], but an often neglected factor is coarser particles 
(PM2.5 to PM10) formed by dust suspensions (from road dust and dust storms) which could also have a 
great influence on human health. 

Nowadays, air quality index (AQI), a numerical model of air quality evaluation based on the pollution 
standards index (PSI) proposed by the United States Environmental Protection Agency (EPA) 
in the 1970s, is widely used to inform the public how polluted the air is. Many countries and regions 
around the world employ different classification methods according to their characteristics. 
An individual air quality index (IAQI) is assigned to the level of each pollutant and the current 
AQI depends on the primary pollutant level (the highest IAQI of the 6 scores); although convenient 
it cannot accurately reflect the particulate pollution state. Many studies based on revised models 
of AQI or fuzzy set theory also face the same situation [7,8]. 

Furthermore, launching a particulate matter concentration monitoring network to determine the 
distribution situation faces phenomenally high costs [7]. The density of the particulate matter 
concentration ground observation network in China is relatively low, and the distribution of 
observation stations is uneven. Remote sensing techniques, which provide consistent measurements at 
broad-scale and frequent time intervals, have been increasingly used to assess surface level of 
particulate matter concentration at high spatial and temporal resolutions [9]. Hence, remote sensing 
retrieval data is introduced in this paper. 

As to population exposure to particulate matters in China, some relevant research has been done in 
specific cities, such as Beijing [7,10,11], Shanghai [12], Wuhan [13], Lanzhou [14], and Hong Kong [15], 
but the spatial differentiation of population exposure to PM2.5 is rarely reported. 

In this paper, we focus on the evaluation of atmospheric particulate matter pollution and propose 
a new particulate air pollution index (PAQI) which only takes PM10 and PM2.5 into account. 
Specifically, we use the retrieved PM2.5 and PM10 concentration over mainland China in 2010 
to calculate the individual indexes of PM2.5 (IPM2.5) and PM10 (IPM10). Furthermore, the spatial 
distribution of PAQI over mainland China in 2010 is discussed. Meanwhile, two indicators of 
population exposure to PM2.5 are calculated and discussed. 

The rest of this paper is organized as follows: in Section 2, we give a brief description of the study area, 
data sources, and methodologies. Then in Section 3.1 we discuss and compare the spatial distributions of 
IPM2.5, IPM10, and PAQI over mainland China in 2010. Section 3.2 discusses Chinese population exposure 
to PM2.5 in 2010. Finally, we summarize our work and draw some conclusions. 
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2. Materials and Methods 

2.1. Study Area and Data Sources 

In this paper, our spatial domain includes mainland China (excluding Hainan Province, which is not 
included in the study area due to the influence of clouds and aerosol type). Three data sources are used 
in the study: 

(1) MODIS atmosphere aerosol product (MOD04L2) is derived from National Aeronautics and 
Space Administration (NASA) at a spatial resolution of 0.1 degree, details on this data can be 
found elsewhere [16,17]. 

(2) Meteorological data (ERA-Interim), the latest global atmospheric reanalysis data produced by the 
European Centre for Medium-Range Weather Forecasts (ECMWF). Dee et al. [18] have covered 
the details of the ERA-Interim product. 

(3) Population density gridded data derived from the Center for International Earth Science 
Information Network (CIESIN) contains the human population distribution. It is a gridded data 
product that renders global population data at the scale and extent required to demonstrate the 
spatial relationship of human population and the environment across the globe. 

2.2. PMjo and PM 2 .5 Concentration Estimation 

Gupta et al. have determined that the artificial neural network (ANN) algorithm can solve the 
PM estimation problem and there are some relevant studies about this topic [9,19,20]. In this paper, 
we make the following improvements to these models: the Levenberg-Marquardt (L-M) back-propagation 
ANN is introduced to build the particulate matter concentration estimation model in this research. 
First, we use stepwise regression to determine the input variables: latitude, longitude, wind speed (WS), 
relative humidity (RH), skin temperature (SKT), boundary layer height (HPBL), aerosol optical 
thickness (AOT), and single scattering albedo (SSA). Among these elements, AOT and SSA are 
obtained from MOD04 L2, and the other four are derived from ERA-Interim. This step can make 
the neural network waste less resources in extracting and fitting irrelevant features. 

The estimation modeling process is described in Figure 1 . We use data of year 2008 to build the 
estimation model and then estimate the PM concentration in 2010. As to the ANN parameters, 
it has been proven that neural networks with a single hidden layer are universal approximators capable 
of representing any real-valued continuous function to arbitrary precision over a finite domain 
if enough hidden nodes are used [21]. However, networks with multiple hidden layers can sometimes 
perform better than single-hidden-layer networks, with fewer total nodes. In this paper, 
the experiment shows that over 16 hidden nodes and 9,000 runs with random sampling provided 
sufficient data samples for the used training data, because the satellite retrieval accuracy of PM 
concentration almost remained the same as the number of hidden nodes and runs increased 
subsequently. Consequently, the ANN estimation model is constructed with eight nodes and 16 hidden 
layers. The accuracy of estimation model will be discussed in Section 3.1. 
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2.3. PAQI Calculation 

Though AQI is calculated on the basis of ground monitoring data, it only reports the air quality 
around each monitoring stations. However, satellite retrievals can exhaustively get the air quality 
spatial distribution, even if it's not accurate enough. Before the PAQI calculation, IPM2.5, 
IPM10 are calculated with the estimation data according to the Chinese AQI standard [22]. 
PAQI calculation is based on a weighting method; weight values are given to IPM2.5 and IPM10, 
respectively, according to a certain condition [23,24], that is, individual indexes and its weight value 
equal PAQI. In this paper, we only consider PM10 and PM2.5 (that is, i = 2), however, it's scalable if 
we have other PM data (PMi, PM0.5, etc.) in the future. First, we calculate the weighting of the z'th 
individual index (Q,). Suppose: 



Q, 



where, 0 < Q, ■< 1 (1 = 0,1,2,...), I, is the z'th individual air quality index. However, if we calculate PAQI 
with Q„ it would be not relate well with the facts under certain circumstances, for example, 
the dust-stormy weather, IPM10 is extremely high while IPM2.5 is relatively low. Hence, we propose the 
weight assignment equation (Equation 2) is: 
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P. = 



0.05 



0 < Q, < 0.05 



LTq, 



(1 - 0.05s - 0.5t) 0.05 < Q, < 0.5 



0.5 



0.5 < Q, < 1 



(2) 



where, s is the number that 0< Q ; < 0.05; t is the number that 0.5 < Q*< 1. P, stands for the modified 
weighting of the z'th individual index. With this step we deal with the extreme values and modify 
the weighting of each individual particulate matter pollutant. Accordingly, PAQI is calculated with 
Equation 3: 



PAQI =5^PA 



(3) 



2.4. Per Capita and Population-Weighted PM 2 .5 Concentration Calculation 

In this paper, per capita PM2.5 concentration (PC-PM2.5) and population-weighted PM2.5 concentration 
(PW-PM2.5) are calculated with Equations 4 and 5 based on grid computing to indicate population 
exposure to PM2.5. Both PC-PM2.5 andPW-PM 2 .5 are obtained from provincial statistics: 

V" PM 

PC-PM 25 = ^ n (4 ) 

2-1 ,=1 i 



PW-PM 25 =Y J (PM i x 



Y" P 



(5) 



where, PM, is defined as the z'th pixel value of PM2.5 concentration, P, is the z'th pixel value of 
population density, n is the total pixel number of the certain province. PC-PM2.5 indicates the average 
exposure level of particulate matter of specific province, while PW-PM2.5 places emphasis on 
the actual effect of particulate matter to residents, taking population distribution into account. 
The weightiness of PM2.5 concentration in a dense population area is larger than in a suburb or sparse 
population area. The PW-PM2.5 method can indicate the PM2.5 exposure level of each resident every day. 

3. Results and Discussion 

3.1. Spatial Distribution of PM Pollution 

Comparing the estimation result with the PM observation data, the comparative analysis consists of 
two parts: accuracies and absolute percentage errors (APE) [19], the correlation coefficients of 
estimated PM10 and PM2.5 are respectively 0.85 and 0.82 (Figure 2), while the absolute percentage 
errors are about 29% and 25%. These results demonstrate the feasibility of the estimation model: 

100x(|y -Y. \\ 
APE = ^ ^ (6) 

Y 0 bs 



where, Y est is the estimated concentration and Y 0 b s is the observed concentration in validation data set. 
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Figure 2. Accuracy of particulate estimation concentrations, (a) The accuracy of PM2.5 
estimation; (b) The accuracy of PMi 0 estimation. 




The spatial distributions of particulate matter pollution over China 2010 are shown in Figure 3. 
The direct comparison between PM10 and PM2.5 concentration is meaningless due to the different 
threshold standards, hence IPM10 and IPM2.5 are calculated according to the Chinese AQI standard. 
Figure 3c and 3d reveal an obvious spatial differentiation. As to IPM10, it presents a good or light 
pollution state, while IPM2.5 shows a moderate pollution situation in the eastern China, 
and a particularly heavy pollution in Hebei Province. Overall, China suffered relatively more serious 
PM2.5 pollution than PM10 pollution in 2010. 

Figure 3. (a) Spatial distribution of PM10 concentration, (b) Spatial distribution of PM2.5 
concentration, (c) Spatial distribution of IAQIpmio concentration, (d) Spatial distribution of 
IAQIpM2.5 concentration. 
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Owing to the spatial differentiation between IPMio and IPM2.5, when evaluating particulate 
pollution level of a certain area or doing research on particulate health effects, both PMi 0 and PM2.5 
should be taken into account at the same time. That's why PAQI is needed. With the PAQI calculation 
method described in Section 3, the spatial distribution of PAQI over mainland China in 2010 is 
calculated (Figure 4). The classification criterion here for PAQI follows the Chinese AQI standard. 

Figure 4. The spatial distribution of calculated PAQI over mainland China in 2010. 
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Compared with the two individual air quality indexes, PAQI shows an entirely different situation. 
PAQI over mainland China in 2010 is overall less than 200 (threshold of moderate pollution). 
There are three regions obviously in a moderate-pollution situation — Hubei Province, 
the Sichuan-Chongqing border region and the Ningxia-Inner Mongolia border region. 

One possible explanation for this phenomenon is that the high value of Hubei Province appeared 
in Jingmen City and its neighboring region may be associated with coal burning and the construction 
dust; the Sichuan Basin is affected by the meteorological conditions, which are not conductive to 
the dilution of atmospheric particulate matter concentration; the high value in the Ningxia-Inner 
Mongolia border region (Shizuishan City, Wuhai City and Dengkou City, China) has a certain relation 
to the development of industry, it suffers from the most serious pollution in the Yellow River Basin 
because many mining industries are grouped on this region. 

3.2. Population Exposure Assessment 

PC-PM2.5 and PW-PM2.5 of China in 2010 are shown in Figure 5. The spatial distribution of 
the results is illustrated in the province-level. The distribution of PC-PM2.5 and PW-PM2.5 is apparently 
different. The high-level of PC-PM2.5 of China in 2010 is concentrated in the Eastern area. 
The provinces with highest PC-PM2.5 are Guangdong, Shanghai, and Tianjin. The 2nd highest places 
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are Shandong, Anhui, and Zhejiang. The high-level areas of PW-PM2.5 are Hebei, Chongqing, 
and Shandong, and followed by Henan, Anhui, Zhejiang, Guangdong, and Guangxi. That is, PC-PM2.5 in 
Guangdong, Shanghai, and Tianjin is high in total, but the potential high pollution area is sparsely 
populated. Although PC-PM2.5 in Hebei and Chongqing is lower than in the high-level area, 
the population in the potential high pollution area is denser. Shandong Province is a special case, 
as it is not only located in the high-level concentration area, but also the population in the potential 
high pollution area is dense. The comparison between two indicators of PM2.5 population exposure 
level indicates the effect of population distribution on the regional PM2.5 exposure statics. 

Figure 5. (a) Spatial distribution of per capita PM2.5 concentration (provincial statistical), 
(b) Spatial distribution of population-weighted PM2.5 concentration (provincial statistical). 
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4. Conclusions 



In this study, PMi 0 and PM2.5 concentrations of mainland China in 2010 are estimated with an ANN 
model, and then IPM. 0 and IPM2.5 are calculated based on grid computing, The spatial distribution of 
IPM10 demonstrates obvious differences with IPM2.5. IPM10 presents a good or light pollution state while 
IPM2.5 shows a moderate-pollution situation in eastern China, with a particularly heavily polluted region 
in Hebei Province. 

On the basis of the two individual indexes, we develop a new particulate air quality index (PAQI) 
based on a weighting method. It is shown that the PAQI over mainland China in 2010 shows overall 
a light pollution state, while some regions suffer a moderate-pollution situation. 

As to the population exposure to PM2.5, the provinces with highest per capita PM2.5 concentration 
are Guangdong, Shanghai, and Tianjin, while the high-level areas of population- weighted 
PM2.5 concentration are Hebei, Chongqing, and Shandong. By comparison, the results demonstrate 
that the population-weighted PM2.5 concentration is more in accord with the actual situation and 
suitable for human health risk assessment. In future research, more attention can be paid to improving 
the spatial resolution of PAQI and discussion of its effects on human health. 
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